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We use de Vaucouleurs’ power-law density-distance relation, to study a hierarchical perturbation 
of the Friedmann universe. We solve the Einstein equation and obtain the density contrast and 
the amplification factor for the perturbation. It is shown that, scale-invariant inhomogeneities 
decay in Einstein-de-Sitter universe. On the contrary, in an open universe, the inhomogeneities 
grow. For low values of S7, amplification peaks sharply and the fluctuations can grow by up to 
a factor of 10 13 from the recombination to the present time. Our analysis of the closed universe 
further confirms that, unlike the common belief that perturbations grow faster with increasing fl, 
scale-invariant perturbations amplify with decreasing It. This result is consistent with the very low 
density expected for a hierarchical universe. 


It is unanimously agreed that, up to a certain scale, 
the distribution of matter in the universe is inhomoge¬ 
neous jlj. However, the scale at which the distribution 
shifts from a power law to a random behaviour is a sub¬ 
ject of strong controversy 11 - For scales up to 100 pc, 
it is well-established that the interstellar medium obeys 
a scaling law jijjT]]. Further up the scale to 150 Mpc, nu¬ 
merous red shift surveys indicate the existence of large 
voids and filaments ||. It is argued that these structures 
have a fractal distribution ||,[| . These facts have brought 
doubts about the hypothesis of the homogeneity, which 
is at the wellspring of the standard cosmology ^ |ll| . 

It remains an ambitious task to study the dynamics 
of a fractal universe. However, analytic inhomogeneous 
universes and also, perturbatively, the evolution of in¬ 
homogeneities in a homogeneous background have been 
studied extensively. De Vaucouleurs was among the first 
to study a hierarchical distribution and the possibility of 
a universal density-distance power law in cosmology [|l 2 | . 
The relativistic inhomogeneous cosmological models and 
specially Tolman spacetime have also been studied widely 
00 - Tolman’s dust solution has been used to model 
a relativistic hierarchical (fractal) cosmology compatible 
with the observational analysis of the redshift surveys 
p5{ . In the perturbative approach, the growth of small 
seeds of inhomogeneities is studied mainly in the con¬ 
text of the structure formation ]I]|i] [H]]. The large scale 
structure of the universe is believed to have grown, due 
to gravitational instability, from small primordial den¬ 
sity perturbations which are fully characterized by the 
density contrast. 

In this article, we attempt to bring these ideas to¬ 
gether. A fractal structure is necessarily scale-invariant, 
eventhough the reverse statement is not true. We con¬ 
sider a scale-invariant spherically symmetric inhomo¬ 
geneous, but isotropic, universe which allows a non¬ 
vanishing pressure. The metric is the Friedmann- 
Robertson-Walker metric with radial-dependent scale 
factor and is contained in the Tolman solution for dust. 


To solve the Einstein equation, we expand the scale fac¬ 
tor around the Friedmann scale factor by a small scale- 
invariant fluctuation. The choice of a scale-invariant 
function for the density perturbation is inspired by de 
Vaucouleurs’ original work on modeling a hierarchical 
universe jl2j]. We take his power-law density-distance 
relation to represent our perturbation and aim at estab¬ 
lishing if such a density fluctuation has any chance of 
growing into large scale structures we observe today. Al¬ 
though the growth of a scale-invariant seed does not im¬ 
ply a fully scale-invariant or a fractal universe, its decay 
would rule out such distributions as viable cosmological 
models. 

To account for the inhomogeneity of space-time, we use 
the spherically symmetric metric, 

diag (g^) = 

(-!, Rp 2 {t,r)r 2 , R P 2 {t, r)r 2 sin 2 d^j , (1) 

where the scale factor R p (t, r) is a function of both time 
and coordinate. The function (1 — fcr 2 )^ 1 multiplying 
R p 2 (t, r) allows for the Friedman-Robertson-Walker met¬ 
ric to be recovered when R p is space-coordinate indepen¬ 
dent. We assume that the scale factor differs from the 
background homogeneous scale factor R(t) by a scale- 
invariant perturbation term and write 

Rp(t,r) =R(t)+SR(t)r~' y , (2) 

where the comoving radial coordinate has been used. The 
value of the fractal codimension 7 has been recently pro¬ 
posed to be about 1 || whereas it was originally believed 
to be approximately 2 [jl!]. Similarly, we expand around 
the homogeneous background matter distribution p with 
the scale-invariant fluctuation to obtain the density func¬ 
tion 

p p {t, r ) = p(t) + 5p{t)r~ J , (3) 

and the radial pressure 


1 



P r (t,r) = wp{t) + aSp(t)r 7 , (4) 

where a is not necessarily equivalent to w. For practical 
purposes, however, it serves to consider an approximately 
adiabatic perturbation where a is almost equal to w. In 
general, a = u> = i> 2 , where v 8 is the velocity of sound 
which is assumed to be unaffected by the perturbation 
10 - 

The expression (g) for the scale factor is chosen because 
it subsequently implies (||) for the density, which has the 
same form as that suggested by de Vaucouleurs Q . The 
radial dependence of the scale factor leads to the same 
radial dependence for the density through Einstein equa¬ 
tion. This scale-invariant behaviour is transfered to the 
pressure through the equation of state. Furthermore, the 
scale factor and the density perturbations are required, 
by the energy-momentum conservation equation, to have 
the same fractal codimension. 

The perturbation of the scale factor, (|j), is a genuine 
perturbation of the metric and not just a gauge mode. 
It is not difficult to see that a solution of Sg^ = 
such that 6go r = 0 and Sg rr ~ t a r~ 7 for a 
general a, cannot be constructed Jl 6 |. 

Our perturbation is only valid for comoving distances 
larger than the correlation length, where the distribution 
is homogeneous and an average density is defined. The 
correlation length has the controversial minimum value 
of 5 Mpc and maximum value of 25 Mpc for some of 
the 3D galaxy catalogues |l] [||[|| . For distances smaller 
than the correlation length, the fluctuation is large and 
the density can be non-analytic with no specific average 
value [||. We are also assuming that crossover to homo¬ 
geneity is reached at some stage. Otherwise, the notion 
of correlation length is meaningless and the present ana¬ 
lytic approach is not applicable ||. 

The Einstein equation is written for the metric (f) 
and the perturbations The zeroth-order contri¬ 

butions, i.e. those with vanishing fractal codimension, to 
the Einstein equation give the Friedmann equations. The 
first-order contributions lead to the second-order differ¬ 
ential equation 

Y + ^- 67 rG(a + 1 )(ap - P) - A 3 Y = 0 (5) 


where 


Y(R) = i?" ( 3 a+ 1 )/ 2 W?, 

A 3 = (1 + 7 ) [1 - a(y - 3)] - ^(a + l)(3a + 1). ( 6 ) 


In the Einstein-de Sitter model, i.e. for k = 0, the dif¬ 
ferential equation (||) has exact solutions which simplify 
to 


SR ~ t x 


3a + l 11 
3(1+ w) + 2 ± 2 


w — 2 a — 1 


u) + 1 


( 7 ) 


at large times. The scale factor fluctuation grows when 
[a + 1) and (u> + 1) have the same sign and 2\a + 1| > 

k + l|- 

The solutions (^) can be substituted back in the time¬ 
time component of the Einstein equation to find the per¬ 
turbation in the density. The density contrast is subse¬ 
quently given by 

m = Sp = t =^ ( 8 ) 

P 

which shows that our perturbation does not grow. Specif¬ 
ically, in matter or radiation -dominated eras the pertur¬ 
bation decays as which is the same as the standard 
decaying mode H [ll[]. Indeed, a simple calculation of 
the amplification factor, 5(to)/6(t), shows that the per¬ 
turbation decays by a factor of 10 -5 in the pressureless 
era (i.e., for t/to ~ 10 5 ). The second solution for the 
density contrast is just a constant. The standard grow¬ 
ing mode of f 2 / 3 in the matter-dominated era and t in 
the radiation-dominated era [ll|] do not arise in our 
scheme. In order to recover these results, we need go 
to higher perturbative orders and expand the scale fac¬ 
tor (and subsequently density and pressure) in series of 
inverse powers of r, i.e. as, 


R p (r,t) 



n =0 


( 9 ) 


where n = 1 term is equivalent to our perturbation for 
7 ~ 1. Had we used this perturbation instead of (2j) 
we would have obtained the standard growing modes at 
the third order of the perturbation (i.e., at n = 3) [^ 8 | . 
Therefore, our analysis shows that, in spite of the grow¬ 
ing modes, the perturbation cannot grow to scales much 
larger than the correlation length, in a flat universe. On 
the contrary, we shall see that in an open universe a grow¬ 
ing mode exists for these scales. 

The differential equation (|s|) for a non-flat universe, 
k ^ 0, is more complicated. To start with, we consider a 
matter-dominated universe where both the homogeneous 
background pressure and the perturbative contributions 
from the inhomogeneities vanish (i.e. a = 0 , u> = 0 ). 
The Friedmann scale factor is given parametrically by 


R = 


d t 
dip 


t = 


0 


2H(1 - H ) 3 / 2 


(ip — sim p), 


( 10 ) 


where ip is real for an open universe and imaginary, ip = 
—id, for a closed universe, and fi is the ratio of the present 
to the critical density 10 - The differential equation 
(^) is now a hypergeometric equation which can be solved 
exactly. The density contrast is given by the growing 
mode 


S+(ip) 


sin (^ 7 ip) 

V7 

+ COS (ypftp) 


—3 sinh '0 
(cosh'i/’ — l ) 2 
5 + cosh ip 
cosh ip — 1 ’ 


27 sinh^ 
cosh ip — 1 

( 11 ) 
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and the decaying mode 


MVO 


cos (V 7 V 0 


Vl sin(v / 7 ^)) 


—3 sinh?/> 
(coshi/’ — l ) 2 
5 + cosh ip 


cosh ip — 1 


2y sinh'!/’ 
cosh^ — 1 

( 12 ) 


We note that for 7 = 0, the above modes coincide with 
the usual results in the literature | 0 . Contrary to the 
decay of fluctuations in the flat universe, scale-invariant 
inhomogeneities can grow in an open matter-dominated 
universe. Infact, unlike the standard open universe den¬ 
sity contrast which remains constant at a redshift of ap¬ 
proximately 2/512 (i.e., cosh?/’ < 5, Fig. |l|) , our density 
contrast oscillates with ip (Fig. |j). 



FIG. 3. The graph of standard amplification factor (7 = 0) 
verses 12 . 




FIG. 1. The graph of standard density contrast (7 — 0) FIG. 4. The graphs of amplification factor at 7 = 1 verses 

verses ip for an open universe. q 



FIG. 2. The graph of density contrast, at 7 = 1 verses ip, 
for an open universe . 

For the amplification factor, 6{to)/d(t r ), from the time 
of recombination to now, once again we have the stan¬ 
dard results for 7 = 0 in which the amplification rises 
monotonically with 12 (Fig. |j). However, for 7 = 1, we 
see that the amplification varies with 12 in an oscillatory 
manner (Fig. |i|). Furthermore, going to a very high res¬ 
olution, at small values of Omega, we see a distinct peak 


at 12 « 0.001 where an amplification by up to a factor 
of 10 13 is obtained (Fig. |J). The amplification is finely 
tuned to the value of 12. For example an amplification 
of 10 5 occurs for 12 « 0.00101629. At higher values of 
12 no peaks and at much lower values of 12 many more 
peaks are observed. This remarkable phenomenon does 
not occur in the standard analysis at 7 = 0. Indeed, a 
self-similar or hierarchical universe is expected to have a 
very low density which is why our perturbation decays 
in a flat universe ||. Here, we have shown that for the 
growth of a scale-invariant inhomogeneity to be in accor¬ 
dance with the anisotropy of the microwave background 
radiation, an upper limit of 12 ss 0.001 is required for the 
density of a hierarchical universe. 

The above results can be verified further by studying 
the density contrast and amplification factor in a closed 
universe. In the standard analysis, the amplification fac¬ 
tor increases with increasing 12 10 - Using equations 
( |lT| ) and ( 0 ) for a closed universe ( i.e. ip = —id), we 
see that the decaying mode ( 0 ) is imaginary for 12 > 1 
and the growing mode 0) now decays. The graph of 
amplification factor shows that the perturbation decays 
significantly from the recombination to the present time 
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FIG. 5. The graphs of amplification factor at 7 = 1 for 
small values of 17. 

(Fig. [6). It also shows that the decay factor increases 
with increasing fl. 



FIG. 6 . The graphs of amplification factor at 7 = 1 for a 
closed universe, 17 > 1 . 


For a matter-dominated universe with a non-vanishing 
perturbation pressure ( i.e. u = 0 but a/0) the differen¬ 
tial equation (||) is no longer a hypergeometric equation. 
However, it can be transformed into a hypergeometric 
equation by using appropriate transformations |l7j. The 
solution is given by the hypergeometric function, 


Y = sin 


X F - 


-3a 


3a + 1 


3a T 1 

a,---a, 


6a + 1 
2 



(13) 


where a = \J\ — 4:A^/2 and 9 = —iip. The second solu¬ 
tion is obtained by replacing a with —1 — a in the above. 
Asymptotically, we obtain 


5R~t x ; x = -|a + ^ \Jl - 4 A 3 , 

S ~ ^-fa-t-ii/ 1- 4A 3 ^ 


For complex values of the square-root term in the den¬ 
sity contrast, above, there is no growing mode since a is 
always positive. For real values of the square root, which 
occur for a > 0.8 or a < —0.4 with 7 = 1, the situation 
is more complex. The terms 1 — 4 A 3 can now be writ¬ 
ten as 9a 2 — 4a — 4 whose square root is always smaller 
than 3a. We, therefore, conclude that asymptotically in¬ 
homogeneities cannot grow in an open universe with a 
vanishing background and a non-vanishing perturbation 
pressure. 

Unlike the preceding cases, an exact solution cannot 
be obtained for a universe with non-vanishing pressure. 
In this case, the Friedmann scale factor is given by the 
expression 

2C 1 i?K 1 +“) (1 3(co +1) 9^ + 5 2 1+3h , 

3(w + l) V 2 ’ 2(3w + 1)’ 2(3w + 1)’ 1 

(16) 

where C\ is a constant. Although, the exact form of 
the above hypergeometric function is not unknown, it 
can be easily shown that it reduces to © for a matter- 
dominated universe and, asymptotically, to a linearly- 
increasing scale factor. In the asymptotic limit, the 
differential equation (§ , for an open universe with 
uj > —1/3, reads 

+ Ct~ 3u - 3 Y + A 3 t~ 2 Y = 0, (17) 

at 1 

where C is a constant. This expression reduces to a con¬ 
fluent hypergeometric equation for F(1/t) after the re¬ 
definition : 

Y = t^ 1±vl - 4A3) eTsdiWF ^ , (18) 

1 2 

where t = C ,1 + 3 “T 1 +3^. In the limit of large times, 
F(1/t) approaches a constant value 0 and we recover 
expression © whose behaviour does not depend on u>. 
For the density contrast, we obtain 

S ~ t-¥+ 3 “+2^ 1 - 4A 3 ) (xgj 

which shows that the perturbation grows for all positive 
values of 3w — 3a/2 or for ui > — 0.3 and positive values of 

a. The growth rate of the perturbation is significant in 
the radiation-dominated era (t/to ~ 10 8 ). In this period, 
the perturbation amplifies by a factor of 10 4 for a « 1/3. 
However, since the baryonic matter cannot grow in this 
era, due to its coupling to the radiation, this result can 
only be relevant for non-baryonic dark matter. 

The main results are summarized in table 1. 


(14) 

(15) 
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equation of state 

P = up + aSpr “ 7 

flat 
k = 0 

open 
k = -1 

matter-dominated 
u> = 0 

decays 

grows 

radiation-dominated 

decays 

grows 

for large times 


Table 1: The table of density contrast for flat and open 
universes in matter and radiation-dominated eras. 

We have evaluated the density contrast for a scale- 
invariant perturbation in flat, open and closed uni¬ 
verses. We have found that, in this specific perturba¬ 
tion scheme, fluctuations decay in the Einstein-de-Sitter 
universe whereas they grow in an open universe. The 
standard results for the growing and decaying modes of 
the density contrast in open and closed universes are con¬ 
tained in our results. The decaying mode for a flat uni¬ 
verse is also reproduced by our method. However, the 
absence of the standard growing mode in the flat universe 
shows that the perturbation grows to larger scales in an 
open universe than in a flat universe. In addition, in our 
analysis, the perturbation does not simply amplify with 
increasing Q. It has an oscillatory behaviour. A remark¬ 
able feature of our perturbation, which is not observed 
in the standard analysis, is the existence of sharp peaks 
at Q < 0.001. At fi « 0.001 the perturbation can be 
amplified by up to a factor of 10 , from the recombina¬ 
tion to the present time. We further confirm our results 
by showing that the perturbation decays even faster in a 
closed universe. We, therefore, conclude that for scale- 
invariant density perturbations the standard belief that 
the density contrast grows faster with increasing f 1 is not 
valid. Our results are fully consistent with the low den¬ 
sity expected for a hierarchical universe. 
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